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Abstract. We show that the particle motion in Bohmian mechanics, given by 
the solution of an ordinary differential equation, exists globally: For a large class 
of potentials the singularities of the velocity field and infinity will not be reached 
in finite time for typical initial values. A substantial part of the analysis is based 

c3 _ 

on the probabilistic significance of the quantum flux. We elucidate the connection 
between the conditions necessary for global existence and the self-adjointness of 



X 



the Schrodinger Hamiltonian. 



1 Introduction 



Bohmian mechanics || f|, |13|, 111, [DJ is a Galilean invariant theory for the mo- 



tion of point particles. Consider a system of N particles with masses m 1; ...,m N 
and potential V = V(Qi, . . . , Qn), where Q fe e IR^ denotes the position of the 
k-th particle. The relevant configuration space is an open subset of uN = d- 
dimensional space H d , for example the complement of the set of singularities of 
V, and shall be denoted by Q. The state of the A r -particle system is given by the 
configuration Q = (Q 1; . . . , Q^r) G Q and the Schrodinger wave function tp on 
configuration space Q. On the subset of Q where the wave function tp ^ and is 
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differentiable, it generates a velocity field = (vf , . . . , v%) 

v fc (?) = — Im — T~r^ — I 1 ) 

m k ip(q) 

the integral curves of which are the trajectories of the particles. Thus the time 
evolution of the state (Qtjipt) is given by a first-order ordinary differential equa- 
tion for the configuration Q t 

^ = «*(<W (2) 
and Schrodinger's equation for the wave function if> t 

where and denote the gradient and the Laplace operator in IR V and the 
potential V is a real-valued function on Q. 

Bohmian mechanics may be regarded as a fundamental nonrelativistic quan- 
tum theory, from which the quantum formalism — operators as observables, the 
uncertainty principle, etc. — emerges as "measurement" formalism. It resolves 
all problems associated with the measurement problem in quantum mechanics 
[[7], ||, f|, [L3|, [14|, [L7|. It accounts for the "collapse" of the wave function, for quan- 
tum randomness as expressed by Born's law p = \if>\ 2 , and familiar (macroscopic) 
reality. For a thorough analysis of the physics entailed by Bohmian mechanics 
see [0, P|, PHI, and [|TJ for a short overview of JEJ. 

Here we are concerned with the mathematical problem of the existence and 
uniqueness of the motion in Bohmian mechanics, i.e., with establishing that for 
given Q and ipo at some "initial" time to (*o = 0), solutions (Qt,ipt) of (0, 
|3|) with Q to = Qo and ip to — exist uniquely and globally in time. (Note 
that Schrodinger's equation (H) is independent of the particle motion, while for 
solving the Equation (^|) for the particle motion we need the wave function ip t .) 
Our first motivation for addressing this problem is the fact that the velocity field 
([!]) exhibits rather obviously possible catastrophic events for the motion: is 
singular at the nodes of if), i.e., at points where if) = 0, so that the solution would 
break down if a node were reached. Furthermore, the solution may cease to exist 
at singularities of the wave function (if it has singularities), at the boundary of Q 
(if it has a boundary), and because of "explosion," that is the escape to infinity of 
a particle in a finite amount of time — events which have analogues in the iV-body 
problem (of gravitational interaction) in Newtonian mechanics. 
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Recall that the problem of the existence of dynamics in Newtonian mechanics 



is notoriously difficult |26| , |12[ . In addition to the possibility of routine collision 
singularities, the iV-body problem with N > 3 yields marvelous scenarios of 
so-called pseudocollisions, where some particles, while oscillating wildly, reach 
infinity in finite time. Examples of such catastrophies have been constructed by 
Mather and McGehee Plfl by Gerver 0, and by Xia [||. While, for the case 
of a "solar system" with small "planetary" masses, Arnold || established global 
existence (and much more) "for the majority of initial conditions for which the 
eccentricities and inclinations of the Kepler ellipses are small," and while Saari 
[j34jl has established global existence for "almost all initial conditions (in the sense 
of Lebesgue measure and Baire category)" for the 4-body problem, for systems of 
more than four particles it is not known whether the initial conditions leading to 
such catastrophies are atypical, i.e., form a set of Lebesgue measure zero — though 
this is certainly expected by most experts to be the case |12| (though not by all 
f25|j). Indeed, apart from obvious scenarios — such as the particles moving apart 
sufficiently rapidly — and those covered by some version of the KAM theorem ||, 
for iV > 5 it cannot, so far as we know, even be precluded on the basis of what 
has so far been proven that this set has full measure! 

It is remarkable that the situation in the corresponding quantum system is 
very different. In orthodox quantum theory the time evolution of the state ip t is 
given by a one-parameter unitary group U t on a Hilbert space 7i. U t is generated 
by a self-adjoint operator H, which on smooth wave functions in 7i = L 2 (Q) is 
given by 

-- 1 2m k 

i.e., Schrodinger's equation is regarded as the "generator equation" for U t . Hence 
the "problem of the existence of dynamics" for Schrodinger's equation is reduced 
to that of showing that the relevant Hamiltonian H (given by the particular choice 
of the potential V) is self-adjoint. This has been done in great generality, inde- 
pendent of the number of particles and for large classes of potentials, including 
singular potentials like the Coulomb potential, which is of primary physical inter- 
est |20], |32fl . We shall discuss the meaning and the status of the self-adjointness 
of the Hamiltonian from the perspective of Bohmian mechanics in Section 4. It 
may be worthwhile to note, however, that the sufficiency of establishing only the 



H = -E k ^—^ + V = H + V, (4) 



however, this example, which is 1-dimensional, involves an infinite number of binary colli- 
sions before the system explodes and thus does not describe a true pseudocollision. 
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self-adjointness of the Hamiltonian for a satisfactory physical interpretation has 
been questioned by Radin and Simon p9 |: "Interestingly enough, while Kato's 
result 'solves' the dynamical existence question in the quantum case, it says 
nothing about the question of x(t) 2 remaining finite in time! From its physical 
interpretation, proof of such regularity property is clearly desirable." 

In Bohmian mechanics we have not only Schrodinger's equation (|3]) to consider, 
but also the differential equation (fj), governing the motion of the particles. Thus 
the question of existence of the dynamics of Bohmian mechanics draws again 
nearer to the situation in Newtonian mechanics, as it depends now on detailed 
regularity properties of the velocity field (H). Local existence and uniqueness 
of Bohmian trajectories is guaranteed if the velocity field is locally Lipschitz 
continuous. We therefore certainly need greater regularity for the wave function 
if) than merely that ip be in L 2 (Q). 

Global existence is more delicate. In addition to the nodes of if), there are 
singularities comparable to those of Newtonian mechanics. Firstly, even for a 
globally smooth velocity field the solution of (0) may explode, i.e., it may reach 
infinity in finite time. Secondly, the boundary points of Q, typically the singular 
points of the potential, are reflected in singular behavior of the wave function at 
such points, giving rise to singularities in the velocity field (|T]).[| 

The problem we address is the following: Suppose that at some arbitrary 
"initial time" (to = 0) the iV-particle configuration lies in the complement of the 
set of nodes and singularities of ipQ. Does the trajectory develop in a finite amount 
of time into a singularity of the velocity field t>^, or does it reach infinity in finite 
time? According to Theorem ^Tl] and Corollary |3.2| , the answer is negative for 
"typical" initial values, for a large class of potentials including the physically 
most interesting case of iV-particle Coulomb interaction with arbitrary charges 
and masses. While we consider in this paper only particles without spin, Bohmian 
mechanics can be naturally defined for spinor-valued wave functions as well j|, 
|8|, [17|. We shall deal with spin, including the motion in a magnetic field, in a 
subsequent work. 

The quantity of central importance for our proof of these results — as well as 
for the question of the self-adjointness of the Hamiltonian — turns out to be the 
quantum flux J^(q, t) = (j^, \ip\ 2 ), a (rf+l)-vector, with = v^\i/j\ 2 the quantum 

2 For example, the ground state wave function of one particle in a Coulomb potential V{q) = 
1/Mj Q £ ("hydrogen atom") has the form e~' 9 ', which is not differentiate at the point 
q = of the potential singularity. 
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probability current. The absolute value of the flux through any surface in f2 x IR 
controls the probability that a trajectory crosses that surface. Surfaces of interest 
for us are the boundaries of neighborhoods around all the singular points for 
Bohmian mechanics. Loosely speaking, the importance of the quantum flux flows 
from the following insight: "If there is no absolute flux into the singular points, 
the singular points are not reached." 

We remark that the quantum flux is, in fact, important for most applications 
of quantum physics, as well as for the mathematics revolving around the self- 
adjointness of Schrodinger operators. Heuristically, the "right" behavior of the 
quantum flux at the critical points ensures self-adjointness of the Hamiltonian — 
i.e., conservation of probability. But suppose we ask, probability of what? The 
usual answer — the probability of finding a particle in a certain region — is jus- 
tified by Bohmian mechanics: A particle is found in a certain region because, 
in fact, it's there. By incorporating the positions of the particles into the theo- 
ry, and thus interpreting the quantum flux as a flux of particles moving along 
trajectories, Bohmian mechanics can be regarded as providing the basis for all 
intuitive reasoning in quantum mechanics. (For more on this point, see also 
[g 1,10, 0, |,|.) 



The paper is organized as follows: In Section 2, the relevant notion of "typical- 
ity" is discussed. Section 3 contains our main results. In Section 3.1 we present 
the broad structure of the argument and in Section 3.2 we show how to trans- 
form the problem to that of controlling flux integrals. The main theorem and 
corollary are proven in Section 3.3. In Section 4 we discuss various aspects of 
the self-adjointness of the Hamiltonian from the point of view of Bohmian me- 
chanics. In particular, in Section 4.1 we show that in d — 1 dimensions global 
existence holds under conditions which in certain respects are milder than those 



of Theorem 3.1 



This is the first work concerned with a rigorous examination of the problem 
of existence of the motion in Bohmian mechanics. For the related theory of 
Nelson, stochastic mechanics, this question has been discussed by Nelson [PB[ 
and also by Carlen M. The behavior of the Bohmian motion at the nodes of ip 



has been addressed by Bohm 0] and Holland | 17fl . Bohm argues that particles 
are either repelled from the nodes or cross them with infinite speed. (Bohm, 
however, was not concerned with the question of existence but with consistency 
with "p = IV'I 2 -") Holland claims to show that a trajectory cannot reach a node 
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unless it is always at some node. His argument, however, is circular, in that it 
requires the very regularity whose breakdown at nodes is the source of difficulty 
Here is a simple counterexample to the claims of Bohm and Holland: Consider 
the one-dimensional harmonic oscillator (with % = m = u = 1) and take as the 
wave function of the particle a superposition of the ground state and the second 
excited state, i/jt(q) = e~ q2 ^ 2 e~ lt ^ 2 [l + (1 — 2q 2 )e~ 2lt }. This wave function has 
nodes (among others) at q = 0, t = (n + \)tt for all integers n. It leads to a 
velocity field which is an odd function of q, i.e., which defines a motion which is 
reflection invariant. Therefore Qt = 0, t ^ (n + |)7r, is a solution of (|2|) which 
runs — first — into the node (0,7r/2) (with velocity and which furthermore may 
be consistently continued through the nodes). 



2 Equivariance and Typicality 

The dynamical system defined by Bohmian mechanics is associated with a natural 
measure, given by the density \4>o\ 2 on configuration space Q. If i(jq is normalized, 
i.e., if the L 2 -norm ||^ || — (In l^o^dq) 1 ^ 2 = 1, then the density \ipo\ 2 defines 
a probability measure on configuration space Q, which we shall denote by P, 
that plays the role usually played by the "equilibrium measure." Thus P de- 



fines our notion of "typicality" []13||. Given the existence of the dynamics Q t 



for configurations — the result we establish here — the notion of typicality is time 



independent by equivariance p3| : 



Po = l^ol 2 Pt = W for all t G H, (5) 

where pt denotes the probability density on configuration space Q at time t— 
the image density of po under the motion Q t . This follows from comparing the 
continuity equation for an ensemble density p t (q) 



dpt(q) 



N 



with the quantum continuity equation 



+ £Vfc-[v?(?)ft(ff)] = (6) 



^ + EVdf(^0 (7) 



and noting that the quantum probability current = (jf, j]v) is given by 

it = ^M 2 = ^-^* v ^ " ^ v ^*)- ( 8 ) 



6 



We further denote the space-time current, the quantum flux, by = (j^, \ip\ 2 )- 
In our proof of global existence, this quantity gives the basic estimate for the 
probability that a trajectory reaches singularities of the velocity field or infinity. 

It is at this stage important to bear in mind the conceptual difference between 
the Equations @ and (0). The continuity equation (|), even without global 
existence of differentiable trajectories Qt, holds "locally" on the set where is 
smooth, with p t suitably interpreted. This understanding is indeed basic to all 
our proofs. 

Equation (|7|), on the other hand, is an identity for every ip t which satisfies 
Schrodinger's equation classically. This is seen by calculating 

^ = ^H^t-AHr t )- (9) 

But, without having established global existence, it is not a continuity equation 
in the classical sense — despite its name. By establishing global existence, we 
simultaneously show that the quantum probability current y is indeed a classical 
probability current, propagating the ensemble density \ip\ 2 along the integral 
curves of the vector field v^. 



3 Global existence and uniqueness 

We make the following general assumptions: 
Al: The potential V is a C°° -function on Q. 

A2: The Hamiltonian H is a self-adjoint extension of H\c°°(n) with domain 
V{H). 

A3: The initial wave function ip is a C°°-vector of H, ip e C°°(H), and is 

normalized, \\ipo\\ — 1- 
The boundary dfl of the configuration space Q will be denoted by S. (Recall 
that usually S is the set of singularities of the potential.) C^(Q), the set of 
C°°-functions with compact support contained in Q, is dense in L 2 (Q), and the 
Hamiltonian is symmetric on this set. Since H is real, i.e., commutes with com- 
plex conjugation, there always exist self-adjoint extensions. The set of admissible 
initial wave functions, C°°(H) = r\^ = iT>{H n ), is dense in L 2 (Q) and invariant 
under the time evolution e ~ ttH / n ^ an( i [ s therefore a core, i.e., a domain of essential 
self-adjointness for 

3 Some special C°°-vectors are eigenfunctions and "wave packets" i/j £ Ran(P[ a ftj), where 
P[ a ,b] denotes the spectral projection of H to the finite energy interval [a, b]. 
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In Lemma |6.1| we show that as a consequence of A1-A3 we may regard ipt = 
e~ ltH / h ipo as being in C°°(f2 x IR) (and thus as a classical solution of Schrodinger's 
equation). Then the velocity field (cf. ([!])) is C°° on the complement of the set 
Af of nodes of if), Af := {(q,t) E Q x IR : if>(q,t) = 0}, i.e., on the set of "good" 
points 

Q := (fi x R) \JV, 

which is an open subset of IR d x IR. Let Gt denote the slice of G at a fixed time t: 
C/t := f2 \ Aft, where Aft '■ = {q E f2 : ^(g, t) = 0}. Then by a standard theorem of 
existence and uniqueness of ordinary differential equations, for all initial values 
(qo, t ) in G there exist r~(q , t ) < t , r(q , t ) > t , and a unique maximal (non- 
extendible) solution Q of (0) on the time interval (r~(g , to), r(q , to)). From 
continuous dependence on initial values, the domain D of the maximal solution 
Q(t] qo,t ), 

D := {(t,q ,t ) : (q ,t ) E G,t E (r~ (q ,t ),T(q ,t ))}, (10) 

is an open subset of ]R d+2 (and Q is locally Lipschitz continuous on D with 
respect to (t, go, to))- Thus r is lower semi-continuous and hence, in particular, 
measurable. Because of the time translation invariance of the theory, we may fix 
to = 0, writing r(go) for r(go,0), with similar notation for r~. Under additional 
conditions on Q and H (see Corollary |3.2| ), we shall show that r(q ) = oo for 
typical go> i- e -> we show that the solution exists globally in time P-almost surely: 

P(r<oo) = 0. (11) 

This is equivalent to 

VT < oo : P(t < T) = 0. (12) 

(Note that by time translation invariance and equivariance, P(r < oo) = for 
all ip E C°°(H) implies that P(r~ > -oo) = for all ip e C°°(H), so that (0) 
indeed implies global existence and uniqueness.) 

3.1 The program 

We view the maximal solution Q t as a stochastic process on Go equipped with the 
probability measure P, i.e., q is distributed according to the probability density 
\ipo\ 2 . The basic criterion for global existence arises from the following properties 
of the maximal solution. The set of limit points L(q ) of the trajectory starting 
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at go (<?* £ i :^ there is a sequence tk,tk — ■> r = r(g ) with lirrifc^ooQ^, = g*) 
is either empty — this is equivalent to linij y T \Qt\ = oo — or nonempty, in which 
case, if r < oo, (g*,r) £ dQ for all g* £ L. (The solution Q need not be 
continuous at t — r, i.e., L might contain more than one point, and there might 
additionally be sequences tk — > t along which \Qt k \ — > oo.) We thus have to see 
whether trajectories come too close to the boundary of Q or to infinity. We do this 
by checking whether they reach the boundary of Q n , an increasing (Q ni C Q n2 
for n\ < n 2 ) sequence of open sets Q n C Q, Q n C Q, which are bounded in 
configuration space, i.e., for all T £ 1R, the set G^t] := G n H (H d x [0, T]) is 
bounded. 

For q £ Qq , we introduce the stopping time r n (q ), at which the process Q t 
first hits the boundary of Q n : 

r n (q ) := sup{s > : (Q t (q Q ),t) £ Q n for all t < s}. 

Now, from the elementary theory of ordinary differential equations, for all go £ 

Qn 
»0 i 

if r(g ) < oo, then r n (g ) < r(g ) and (Q T ^ qo )(qo), r n (q )) € dQ n . (13) 

Furthermore, the sequence r n is increasing in n. For all n and all T < oo, we 
have 

{q eg :r<T}c (G \ G%) U {g G ^" = ^ < T} 

and therefore 

P({?o eG :r<T})< P(G \ g%) + P({g £ Ql : r" < T}) (14) 

Thus to obtain the global existence and uniqueness of Bohmian mechanics for 
typical initial configurations, it is sufficient to establish the vanishing of the right 
hand side of fli~4]) as n — > oo for some sequence of sets C/ n . (Note as a matter of 
fact that the right hand side of (13) decreases as n increases.) 



To proceed, we need to separate different parts of the boundary of Q n which 
we shall treat in different ways: those close to infinity, those close to S — dQ, 
and those close to the set Af of nodes of the wave function. We introduce /C n , a 
sequence of bounded open sets exhausting JR d , K, n / !R d ; S 6 , a sequence of closed 
neighborhoods of S of "thickness 5" ; and Af e , a sequence of closed neighborhoods 
of Af of "thickness e." (For the following general remarks, we do not need to 
specify these sets more concretely; this will be done in Section EO|) Thus the 
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index n accordingly gets replaced by e8n. Q e5n then denotes the set of "e-5-n- 
good" points in configuration-space-time: 

g* Sn : = (((/c n nn)\ s s ) x ir) \ j\f e , 

and Qf n denotes the slice at a fixed time t 6 H: 

ge6n._ (/C"nfi)\(5 5 UAT/). 

Furthermore, we define 

$ T) :=((rnQ)\5 5 )x(0,T). (15) 
From (0), we may write, with x := (Q m i n (T eSn ,T), min(r e<5ri , T)), 

e £o 5n : ^ £<5n <T} = {g e : x e dg e5n n (IR d x (0, T))} 

u {g e $f n : x e {dS 5 n ft) x (o, T)} 

u {go e gf n : x e (<9/C™ n fi) x (o, T)}. (16) 

and therefore we arrive at 

P({?o e : r e5n < T}) < P({g e g e 5n :xedM e n g$ tT) }) 
+ P({?o e £o 5n = * e (as 4 n ft) x (o,T)}) 
+ P({?o e : x e (dlC n n ft) x (o,T)}). (17) 

By virtue of (almost sure) global existence follows if for some suitable choice 
of sets J\f e , S s , and fC n , P(Go \ G(f n ) and the right hand side of ( [P7D can be made 
arbitrarily small by appropriately choosing e, 5 and n. 

3.2 The flux argument 

Consider the random trajectory (Go,P,Qt) obtained by stopping the original 
process Q t at time r and placing it in the cemetery f : The process Q t : Go — ► 
ft U {f} is defined, for all t > 0, by 



Qt(g ) for £ < r(g ) 
f for t > r(g ). 



<9t(?o) := 

Let p t be the image density of Q t restricted to ft. 
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Denote by X the set 

1- ■= {{Qt{qa),t) : t E (r~(g ),r(g )) and q £ Q }, 

and by X t := RanQ* \ {f} (It C its time-t slice. X is an open subset of Q. (X 
can be identified with D n ({0} x IR d+1 ), cf. (0).) Clearly ft = on Q t \X t for 
i > 0. Note that on X both |^ t | 2 and p 4 are solutions of the continuity equation (|6|) 
restricted to X with the same initial data. Uniqueness of solutions of quasilinear 
first order partial differential equations on the set where the characteristics exist 
implies that for all t > 

Pt(q) = \Mq)\ 2 for aUgeX t . (19) 

Consider now a smooth surface £ in Recalling the probabilistic meaning 
of the flux Jt(q) '■= (pt(o) v (o) > Pt(o}) > we obtain that the expected number of 
crossings of £ by the random trajectory Q t (including tangential "crossings" in 
which the the trajectory remains on the same side of E) is given by 

/ \J t (q)-U\da (20) 

where U denotes the local unit normal vector at (q, t). (J s ( J-U)da is the expected 
number of signed crossings.) (Consider first a small surface element which the 
trajectory can cross at most once. The probability density for this crossing is 
readily calculated to be | J • U\. Invoking the linearity of the expectation value 
yields then the general statement. )Q 

The probability of crossing S (at least once) is hence bounded by (^(]). From 
(|19D we obtain that 

\Jt ■ u\ < |(|^ t (g)|V%), W) ■ u\ = |(/s |^l 2 ) • u\ = I J* • u\ 

and thus we arrive at the bound 

P(Q t crosses E) < j \J^(q) ■ U\da. (21) 

If now the sets W e , S 5 , and JC n are choosen in such a way that their boundaries 
are piecewise integrable surfaces, the events on the r.h.s. of (|TT| ) are crossings by 

4 In stochastic mechanics |28), which involves the same quantum flux, the particle trajectories 
are realizations of a diffusion process and are hence not differentiable, i.e., velocities do not 
exist. Thus in stochastic mechanics the flux does not have the same probabilistic significance 
and hence the subsequent arguments are not valid for stochastic mechanics. 
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Qt through the respective surfaces, and hence fl2"T|) implies the following bounds 
for the terms in (ITT 



P(^(^ £ n^ T) )) < / \J+{q)'U\da:=N{e,6,n), 

P(xe((dS 5 nQ)x(0,T))) < [ \J*>(q)-U\da:=S(S), (22) 

J(ds s nn)x(o,T) 

P(x G ((dKT 1 flfi) X (0, T))) < I \J^{q) -U\da =:l{n). 

J{dK n n£l)x{0,T) 

(If a boundary happens to be the empty set, the corresponding integral of course 
vanishes.) 

It seems intuitively rather clearf] that all the flux integrals should vanish in 
the limit e — > 0, 5 —>■ 0, and n — > oo: It seems fairly obvious that the "nodal 
integral" N(e, 5, n) should vanish as e — > since J^* is zero at the nodes.0 The 
"singularity integral" S(5) should vanish in the limit <5 — > if the set S has 
codimension greater than 1, which is usually the case. Furthermore, j = at 
S is a natural boundary condition defining a domain of self-adjointness of the 
Hamiltonian. Finally, the "infinity integral" I(n) should tend to zero as n — > oo 
since ipt(q) (which is sufficiently smooth) and hence J^'(?) should rapidly go to 
zero as |g| — > oo. 



3.3 Global existence of Bohmian mechanics 

Our main result is the following theorem: 
Theorem 3.1 Assume A2, A3, and further 

AT: Al and S C UHLi^z? where m < oo and the Si are (d — 3) -dimensional 

hyperplanes; 
A4: J T ||V^t|| 2 ^ < oo for all < T < oo. 
Then P(r < oo) = 0. 

Since <S; is a (d— 3)-dimensional hyperplane, it may be written as Si = {yi = a^} 
with yi denoting the map H d — > 1R 3 , q i— >■ (g • y^ 1 , g • y z 2 , g • ?/f ) where ?/; 2 , yf are 
3 orthogonal unit vectors normal to the hyperplane Si and a; G H 3 a constant. 

5 By mentioning these heuristics we do not wish to suggest the structure of the rigorous proof 
given in the next section, nor need this proof sustain these heuristics. 

6 One might worry about the "size of dj\f" being uncontrollably large. However, since ip is a 
complex smooth function, Af might be expected to have codimension 2 "generically," so dN e 
should have small area. 
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The Condition Al' on the shape of S fits well with the 3-dimensionality of 
physical space. If V is a central potential, <S; is of the form {qj = 0}, and for a 
pair potential, <S; is of the form {q^ — q^- = 0} for some 1 < % < j < N. (Note 
that if d = vN < 3, Assumption Al' demands that S = 0.) 

Under the Assumption Al', the configuration space f2 = 1R q! \S; in particular, 
L 2 (Vl) = L 2 (lR a! ).[] Recall that Hq denotes the self-adjoint operator 

N h 2 

fc =i 2m fe 

on the Hilbert space H = L 2 (VL) = L 2 (lR d ). 

The Condition A4 of "finite integrated kinetic energy" may be ensured 
by bounding the quadratic form (Vipt, V^i) < M(ip t , Ho4>t) with M = 
(2/h 2 ) max(mi, . . . ,mjy) by the form (i/j t , Hiftt), which is finite and independent 
of t for ip (and hence ip t ) in the form domain []3T| Q(H)(d T>(H)) of the Hamil- 
tonian The following corollary shows that Theorem |XT] indeed implies the 
global existence and uniqueness of Bohmian mechanics for all ipo G C°°(H) for a 
large class of Hamiltonians. 

Corollary 3.2 Assume 

Al": Al' and V = Vi+V 2 , where V\ is bounded below, and V 2 is H -form bounded 
with relative bound a < 1, 



A2': H is the form sum H + V 

and A3. Then P(r < oo) = and Bohmian mechanics exists uniquely and 
globally in time P-almost surely. 

Proof. We show that A4 holds: That V% is E.§-ioxm. bounded means that 
Q(H ) C Q{V2) and that for ip E Q(H ) there exist constants a, b > such 
that 



7 Thus Theorem 3.1 does not cover the case of a bounded configuration space f2, for which 



boundary conditions of Dirichlet or Neumann (or mixed) type are normally imposed. See, 



however, our Theorem 4.1. 

8 Note that the notation (?/>, Atp) for the quadratic form associated with the self-adjoint 
operator A is symbolic: Only for ip S T>(A) does it coincides with the indicated scalar product 
in Ti = L 2 (TR d ); more generally it can be defined via the spectral representation for A. 
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Since Vi(q) > -c, c> 0, for all q G Q, we obtain for ^ G Q(#) = Q(if Q ) n Q(Vi) 
that 

(l-o)(V,W) < (V>,(#o + W) + &(^) 

< (ij,(H + V 1 + V 2 )ij) + c(ij,ij) + b(ij,ij) 

= (iP,Hil;) + {b + c)(i>,i>). 
Hence with a < 1 we have that for ip G Q(-ff) C Q(Hq) and all £ 

M 1 — a 1 — a 

' .(^ ,^o) + ^||^o|| 2 



1 — a 1 — a 

and A4 follows. □ 

The class of H -foTm bounded potentials, with arbitrary small relative bound a, 
includes for example R+L°° or L 3 / 2 +L°° on IR 3 , where R is the Rollnik class. (For 
details, see for example |^T], |36|, |52]|.) Therefore such if -f° rm bounded potentials 



include power law interactions l/r a with a < 2, and thus the physically most 
relevant potential of iV-particle Coulomb interaction with arbitrary charges and 
masses. (The class of H -foTm bounded potentials contains the more familiar 
class of iJo(-°P era t° r ) bounded potentials, which already includes the iV-particle 
Coulomb interaction |p0[ .) Furthermore, harmonic and anharmonic (positive) 
potentials are included, and arbitrarily strong positive repulsive potentials. 

Proof of Theorem We establish (O— for all < T < oo, P(r < T) = 0— 



following the program described in Section 3.1 and the flux argument of Section 



3.2 



We first choose suitable sets J\f e , S 5 , and fC n . Let e > 0. Set 

Af*:= |J C*(k), (23) 

where (C e (fc))fceiN is a "partition" of configuration-space-time into closed cubes 
with side length e whose edges are parallel to the canonical basis vectors of 
]R d+1 . Let 5 = (6i,...,8 m ), S t > for all I. Recalling that S C (J™i«S; with 
5/ = {y* = a z})> set 

m 

<S 5 := |J Sf 1 , Sf 1 := {g G IR d : dis%<S z ) < 5,} = {|y, - a,| < 5$. 
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For the cutoff at infinity we choose open balls with radii n G 1R : 

KJ l := {qeJR d : \q\ < n}. 
By virtue of 0, (0), and fl2|), we obtain that for all < T < oo 
P(r<T) < P(Q \Q e Sn ) + P(T e5n <T) 

< P(Go \ Gf l ) + P(* e n £$ T) )) 

+P(a; G x (0, T))) + P(ac G ((<9/C n n Q) x (0, T))) 

< P(eo\0 + N(e,^) + S(5) + I( n ) (24) 

For the first term on the right hand side of (f24"D recall that Q§ n = (/C n D fi) \ 
(A/" e U S 5 ); therefore 

Go \ Gt = (Go \ /C n ) u n S 5 ) u (£ n A/" e ), 

and thus 

P(£o \ Gt) < P(G \ K- n ) + P(G n S s ) + P(Go n A/^). 

The vanishing of the three terms on the right hand side in the limit n — > oo, 
8 — > 0, resp. e — > 0, follows easily from the facts that P is a probability measure 
with density \ipo\ 2 , and that the respective sets tend to P-measure sets. 

The vanishing of the remaining terms in fl2"4|) is the content of the following 
lemmas: 

For all < T < oo there exists a sequence rik, 
lim I(rijfc) = 0. 

Lemma 3.4 Assume AV, and A2-A4. Tnen tnere exists a sequence of m-vectors 
_> o as jfe -> oo (usto 3 > /or all l,k), with 

lim S(o" (fc) ) = 0. 

Lemma 3.5 Assume A1-A3. For all < T < oo, n < oo and <5 > 0, 

limN(e,5,n) = 0. 

e^0 y ' 



Lemma 3.3 Assume A1-A4. 
nk — > oo as k — > oo, with 
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These lemmas will be proven below. Lemmas 573], |3.4|, and |3.5| imply that the 



r.h.s. of (24) can be made arbitrarily small. (Note that if d < 3, Assumption Al' 
demands that S = 0, and hence that S = 0, so that Lemma |3.4j is trivial in this 
case.) □ 



Proof of Lemma [57 

I(n) = / I J^Hq) ■ U\ da = [ T ( \f t (q)-u\dsdt 

J{dK n nn)x(o,T) Jo JdK n nn 

< fi [ T I |VV>t| dsdt =: fA{n) 

Jo JdK n nn 

with fi = h/ min(mi, . . . , uin), ds the (d — l)-dimensional surface element of 
dJC n , and u the local unit normal vector of this surface. To show that I(n) goes 
to along some sequence nk, we prove a stronger statement, namely that I(n) is 

poo _ 

integrable over n. This is immediate since / I(n) dn yields the volume integral 

Jo 

of \tp t \ |VV>t|, which is easily estimated: 



L 



l(n)dn = f f \i/) t \ \Vip t \ dqdt 

Jo Jn 

< / ||^t|| l|V^|| * = / ||V^||*<oo, 
Jo Jo 



where we have used A4 for the last inequality. We may thus conclude that there 
exists a sequence {rik)k with rik — > oo as k — > oo, along which I(njfe) — > 0. This 



proves Lemma |3.3| . □ 



Proof of Lemma \3.4 - We may assume that d > 3. We shall use the following 
Inequality: For tp G Q(H ) 

I d JuT^dq< f d m\ 2 dq. (25) 

This is a straightforward extension of the inequality known as Hardy's inequality 
or the "uncertainty principle lemma" (see, for example, |32]|) usually given for 
i) G C^IR 3 ): 

(One immediately obtains ( p5| ) for d = 3 and ^ G C^°(IR 3 ). Then, viewing ip G 
C^ > °(]R' i ) as ^ G C^°(1R 3 ) by keeping all coordinates fixed except y/, one extends 
this inequality easily to C™(1R d ). It is then further extendible to ip G Q(H ) 
because C^°(lR d ) is dense in Q(Hq) with respect to the i?o-form norm.) 



16 



First we estimate 

r , rT r 

ds dt 



S{6) = f \J^(q)-U\da = f f 

JdSSx(0,T) 1 v ' 1 Jo Jes* 

< / / 4l 1^1 |V^ t | dsdt =: /i£S^)- 



i=1 Jo Jas/nn l=1 

We now integrate (l/|yz — aj|)Sj(5{) over 5i = \yi — a z |: By the definition of 
<Sf ! = {|y; — a; | < 5i}, this yields the volume integral of (l^l/lyz — |VV>|, 
which may be bounded as follows: 



Jo Jo Jn yj - a z | 



= / /, i i lWt| dqdt < f || , - H HIV^II^ 

Jo Ju d \yi - a z | Jo |y/ - a z | 

< 2 f \\Vipt\\ 2 dt < oo 
Jo 



using Schwarz's inequality and the Inequality (^). Since 1/5/ is not integrable 
at 5i = 0, for each I there exists a sequence 5^ with 5^ — > as k — > oo, along 



which Sj((J, (fc } ) -> 0. This proves Lemma |3.4| . □ 



Proof of Lemma \3. 5| . This proof is more involved than the previous ones, since 

the nodal set is unknown. The basic idea is the following: Where the (d + 1)- 

dip 

gradient ip' = (Vt/>, ls sma U the current is very small, and where ip' is not 
small the surface area can be controlled. 

Let rj > 0. We split the part of N e contributing to the surface dJ\f e fl £/{ut) 
into two (not necessarily disjoint) sets: 

My := |J C £ (A;), and 7V< := (J C e (A;) with 
feG/> fce/< 

/> := {fc : C e (fc) n {(q,t) : ^( g ,t) = 0, |^(g,f)| > 7?} n G$ T) ? 0} 

/< := {k : C e (A;) n {(g,t) : ij>(q,t) = 0, |^'(g,t)| < V } n £$ T) ^ 0} 

Then 

N(e,5,n)=[ k \J^(q)-U\da < [ \J**(q)\d*+ f \J^(q)\ da (26) 



On the compact set ^S 2 ^"^ (cf. fll5l)) there exist a global Lip schitz constant 



L for ip', and a global bound K for |^)'|. Observe that for e < min(<5/ (2\/rf), l/\/d), 

,(5/2)(n+i; 
'(-1,T+1) 



A/5 C G^li?. Let therefore e < min(5/(2 v / rf), 1/W). 
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Consider first 7V< . In this set the flux \J^\ is very small. We may estimate the 
integral by simply taking an appropriate bound of | J^\ times the total area of the 
surfaces of all the cubes. In every e-cube C e of jV< there is a point (q*,t*) e J\f 
with \i/j'(q*, t*)\ < T]. Thus (in every e-cube of N< and hence) for all (q, t) G 7V< 



W(q,t)\ <v + LVd+le. 

\ip\ is thus bounded on (every e-cube of jV< and hence on) A/"< by (77 + 
+ Ie)V" + le = : Ci//e + c 2 e 2 . The flux is then bounded by 



1^1 = Vd^iT + bi 2 < \ip\ 2 +\f\<\ip\ 2 +n\ip\m (27) 

< (cir/e + c 2 e 2 ) 2 + /i(c^e + c 2 e 2 ) (7; + Wd + le) 



To bound the surface area of A/*<, we simply add the areas of the surfaces of all 

,(«V2)(n+i; 
'(-l.T+l) 



e-cubes in ^/l^™^. The number of e-cubes in GfliT^y is bounded by c 3 /e d+1 



with 

c 3 (n,T,d) = (T + 2)(2n + 2)\ 

and the surface area of a single cube is equal to 2(d + l)e d . Thus for the surface 
area of cW< we have the bound 

|^| < (28) 



and combining Q27D and (28) we obtain that 



|J^a< suplJ^I (IQA^I) 

V^< / 



< 2(rf+1)C3 ((c 1 r ? e + c 2 e 2 ) 2 + /i^e + c 2 e 2 )(7? + Lv'd+Te)) (29) 



0( V 2 ), e^O. 



Consider next the set 7V>. On this set we can control the size of the nodal 
surface. To do this we use a further partition of configuration-space-time into 
cubes (C 7 (fc))fc e iN of side length 7 (with sides parallel to the sides of the C e -cubes). 
We choose 7 so small that any C 7 -cube which contains or overlaps the interior of a 
C e -cube of J\T> lies completely in G^it+i)^ > ^ e -' 7 < mm (^/(2v^)— e, (1/ Vd) — e). 
(7 will later be chosen to be proportional to rj, and we shall take the limit e — > 
for fixed 77, so that 7 ^> e eventually.) 
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We show now that in each 7-cube the number of e-cubes in 7V> is small, at 
least compared with Consider a 7-cube C 7 (A;) containing or overlapping 

the interior of a C e -cube of A/*>. Then there is a point (q*,t*) G Af D C](k) with 
> ^, where C^ik) is the "e- fattened" 7-cube, i.e., the cube of side 
7 + 2e with the same center as C 7 (fc). That ^'(g*, > rj implies that 

for either i = 1 or i = 2 (or both), with ^ := Re?/> and ip2 '■= Im^. 

Let efc be that basis vector which is closest to the direction of ^ (<?*,£*), i.e., 
for which \e k • ^(f?*, is maximal. Thus 

/2(d + T) 

and hence we have that for all (q,t) G C](k) 

V 



\e k -^(q*,t*)\ > 



e k -4(q,t)\ > , ( -L v / rf+l(7 + 2e). 
2(d+l) 



Now choose 7 such that Lyd+1(7 + 2e) = 77/ (2y2(d + 1)J, i.e., introduce 
c 4 := 1/ (2L v / 2(rf + 1)) and set 7 = c 4 t? - 2e. Then for all (q, t) G C 7 (fc) 

\e k -i// i (q,t)\> , 77 . (30) 



2^2(^+1) 

Let a; and y be two space-time points in C]{k) H A/" e with = and 

x — y = le k ,l > 0. Then, on the one hand, by the global bound K on |^'| we 
have that 



\ipi(x)\ < KVd + le. 

On the other hand, it follows from (^TJ) that \ipi(x)\ > lrj/(2^2(d + 1)). Thus 
/ < 2i^v / 2(^ + ^e/ 7 ? = : c^e/r]. Therefore the number of e-cubes in A/"> contained 
in C](k) and lying in an e^-column — the set of e^-translates of an e-cube — is 
bounded by (05/77) + 1. (This is a rather crude estimate. The number of such 
cubes is in fact bounded by 2d + yd + 2, independent of 77, as can easily be seen 
by controlling also the projection of ip[ orthogonal to e k .) 

Now the number of e^-columns in C](k) is no greater than [(7/e) + 2] d , while 
the number of 7-cubes in Q^xx+i) * s bounded by c^/^i d+l . Thus we obtain a 
bound for the surface area of A/">: 
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\J^\ may be estimated (as in (|2?D) by invoking now the global bound K for 
l^'l to yield 

\J^\ < K 2 {d+ l)e 2 + fiK 2 Vd+le 
on Thus we arrive at the estimate 

/ \J\da (31) 

- (f + 1 ) (t)" ^ -V» 2(d + 1>e '( /f 2(t * + 1>t2 + 2 ^ TTt ) 

Using fl29|) and (|3T]) , by letting first e — > and then 77 — 0, it follows from (f26|) 
that limN(e,5,n) = 0. □ 

3.4 Remarks 

3.4.1. It is an immediate consequence of continuous dependence on initial con- 
ditions for solutions of ODE's that the probabilistic negligibility of the set of 
"bad" initial values B := {q G Go : r(g ) < 00}, P(£>) = 0, implies the negligibil- 
ity of B in the topological sense: B is of first category in Q , i.e., it is contained in 
a countable union of nowhere dense (in Q ) sets. (Take B t = {q G Go : r(g) < t}; 
cf. also [0.) In other words: Global existence of Bohmian mechanics is typical 
and generic. 

3.4.2. Since P is equivalent to the Lebesgue measure L on Q , we have also 
that L(£>) = and we thus have the global existence and uniqueness of Bohmian 
mechanics L-a.s. on Qq. 

3.4.3. The flux argument shows that any given hypersurface in Q x 1R (where 
ip is C°°) of codimension greater than 1 will (almost surely) not be reached. 



3.4.4. We have shown that under certain conditions on the initial wave func- 
tion and the Hamiltonian, particle trajectories exist as solutions of (0) globally 
in time for P-almost all initial conditions. In the introduction we have already 
given an example showing that in general (i.e., assuming merely the conditions of 



Theorem 3A or Corollary |3.2|) this result does not hold for all initial configura- 
tions. However, in that example the dynamics is uniquely extendible to a global 
dynamics Q : H 2 — > H, (q, t) 1— > Qt(q)- There are 3 continuous trajectories which 
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periodically run into nodes of the wave function, while the other trajectories are 
global solutions of (H). This extended dynamics Qt(q) is continuous. 

However, if the trajectory running through the node at t = 0, q = 1 is analyzed, 
one finds that locally Qt(l) ~ \J\^ 2 + 1> i- e -> the map Qt(g) is not different iable 
with respect to t at t — for fixed q — 1. This may, for example, be seen by 
considering the flux through q — 1 for t near 0, or, what amounts to the same 
thing, by employing the Formula fl4l|) (see Section ^) expressing the trajectories 
as curves of constant value of the function F(q,t) = \ip t \ dx. (This behavior 

J — oo 

of trajectories hitting nodes is in fact typical — though it does not occur in the 
example for the trajectory at the origin; in fact, if ip(q,t*) has a node of order 
k at q*, ijj(q,t*) ~ ax k with x = q — q*, then F(q,t*) ~ F(q*,t*) + ax 2k+1 , a = 

Jf^, and -fo(q,t*) = -J A *(q) ~ ^ so that F{q,t) ~ F(?*,t*) + ax 2k+1 + 
6x 2fc s + cs 2 , s = t — t*, in the vicinity of the node. Thus for c ^ 0, the equation 



F(q,t) = F(q*,t*) implies that x ~ 2fe4 \/ _ a s2 -) 

Concerning the regularity of Qt(q) in q at fixed t, one sees in the example that 
for suitable choices of initial time the solution map will fail to be different iable 
at q = (where there will be a fifth root singularity) or at q = ±1 (where there 
will be a cube root singularity) as a function of q for fixed t. 

For an even stronger breakdown of regularity in q for fixed t, consider the 
harmonic oscillator in 3 dimensions, and take the (n = 1,1 = l)-state ip{q,t) = 
re -(r 2 +z 2 )/2 e ui> e -5it/2 j n C yii nc i r i ca i coordinates. This wave function vanishes only 

at r = 0, i.e., on the z-axis. Particles circle around the z-axis with angular 
velocity 1/r 2 . The map Q is uniquely extendable to a global dynamics given by a 
continuous map, which is however not different iable with respect to q, by defining 
Qt(qo) = qo for all t and q 6 Mo- 
lt is possible also to give an example in which the extended map must fail 
even to be continuous with respect to q for fixed t: Consider free motion in 
1 dimension, and let the wave function if> be even, (real and positive), C°°, and 
supported on [—b, —a] U [a, b] with < a < b < oo. Then ip e C°°(Ho). Moreover, 
there is a h > such that (e iHotl / h ip^j (0) ^ 0. Let ip = e iH ° tl/h ip. ij} t is then 
even for all t, so that the velocity field is odd, i.e., symmetric under reflection. 
Any extension Q which respects this symmetry must have Qt(0) = for all t. 
Then the map Q is discontinuous in q for t = t\, and, in fact, any extension must 
have this discontinuity. 
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3.4.5. It is well known — at least if V is real analytic in f2 (see for example [^C 
page 98) — that if ip vanishes on a nonempty (bounded) open set in configuration- 
space-time, it vanishes everywhere (in the components of Q x 1R that intersect this 
set). We remark that under the hypotheses of Corollary [372], the same conclusion 
would in fact obtain merely if ip were to vanish everywhere on the boundary of 
such a set (and even with the possible exception of a single piece of the boundary 
contained in a constant-time hyperplane), since it would then follow from global 
existence and the inaccessibility of the nodes that ip must vanish everywhere in 
this set. 

3.4.6. The probability of reaching the nodes P(x G (dAf e fl G?q T \)) may also be 
estimated without using flux integrals. We include this argument, which involves 
a choice for M e different from the one used earlier. We remark that for the new 
M e we can see no reason why dM e must be smooth, even piece-wise. Notice also 



that Lemma 3.G involves both stronger premises and, since the convergence in it 



is uniform, a stronger conclusion than the corresponding Lemma p~5 . 
Lemma 3.6 Assume A1-A4 and, for e > 0, let 

Af e :={(q,t)enxlR: \ij> t (q)\ < e}. (32) 
Then, uniformly in 5 and n, with x := (Q min ( T ^n T ), min(r <E<5ri , T)) 7 
limP({g e Gt : x G {dN e fl £? n T) )}) = 0. 



The proof involves a fairly standard "existence of dynamics" argument and is 
analogous to that of Nelson |28] for the similar problem in stochastic mechan- 



ics: One looks for an "energy" function on the state space of the motion which 
becomes infinite on the catastrophic event. With good a priori bounds on the 
expectation value of that function, one can control the probability of catastrophic 
events. 

Proof. The function which recommends itself here is log i.e., what we control 
is the "entropy." 

We first present a formal estimate, disregarding the problem that the solution 
curve Qt{q) starting at q may not exist for all times — which is taken care of below. 
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Let E denote the expectation with respect to P. We compute for arbitrary T: 

d 



E (| log |Vr(Qr)| - log |Vo| I) = E 

|2 



o dt 



\ogMQt)\dt 



E 



T (l 1 d\MQt)\ 2 | WjHgj ^roVU 
o 1 2 |^(Q.)l a & ltM<WI WtJ J 



< 



E - 



1 1 



o \2MQtW 



d\MQt)\ l 



dt 



dt + J^E Ll 



\vMQt 



MQtW 



dt, 



(33) 



where we used for the inequality the bounds 

V|V>| < |VV>| and \v^\ <fi 



(34) 



Now use the equivariance of |-0| 2 (cf. ©) to compute the expectation E(f t (Q t )) 
In \'4 , t{ ( i)\ 2 {ft{ ( i))d(l and obtain that the right hand side of (|33| ) is equal to 

ami)? 



71 

o Jn 2 



dt 



dqdt + fi f f \Vip t {q)\ 2 dq dt. 
Jo Jn 



(35) 



By virtue of (^) we replace \d\ip t {q)\ 2 /dt\ by \ip^(q)Hip t (q) — i/j t (q)Hi/)%(q)\/h. By 
Schwarz's inequality, the first term of (|35"D is then bounded by 



I- [ T \\M\\HM dt = | ||FVo|| < oo, 
n Jo n 



and the second term is bounded for each T < oo by Assumption A4. 

To construct from this a rigorous proof we need only define a suitable killed 
process. For t > we define Qf n : Qf n U {f} — ► Qf n U {f} by 



Qf n (q) 



Qt{q) for t<r eSn {q) 
f for t > T e5n (q) 



(36) 



For completeness, we set Qf n {\) = f for all t > 0. Consider the probability 
measure Pq* 5 ™ on ^q* 5 ™ U {f } which has the density 

Pt(q) := \Mq)\ 2 ior q e Gf n 

(and, of course, Po Sn (\) = 1 — fgedn p^ n (q) dq). The image measure of the process 
Qf n is denoted by Pf n := Pf n o (Qf and has the density pf n on Qf n . From 
the definition of AT 6 (B2|), 



{go e Sf 1 = x e (dAf n Sfo>))} C {g e ^ n = 1^)1 = e}- (37) 
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Since we keep 8 and n fixed, and since the estimates are independent of 8 and n, 
we will omit the indices 8 and n on Q eSn , Q eSn , p eSn . 
Define for 96^ and t > 



D e t {q) := \og\ip min(T e (q)tt) (Q min(T c {q):t) (q))\ - log |^o(?)|- 



One has that 



D e M = [ J^fGz) dt = f t o Qt(q) dt 



where 







ft(y) ■■- 



for y = I 



1 1 



lw)l 



2|^(y)| 2 0t 
We shall show that uniformly in e 

P({q G ^ : \D e T (q)\ > K}) -> as AT -> 00 



(3« 



Then, since for g as in (ft7|) D^(q ) = loge — log (^((Zo)!, the lemma follows from 
(|38|) by observing that 

P({geS e : |log|Vo(g)||>^})^Oas^^oo 

holds uniformly in e, which is immediate since the density of P is l^o) 2 - 
By Markov's inequality we obtain that 

P({qeF :\D* T (q)\>K})<±EUge j*f t oQldt\. (39) 

Recall now that P = Pq on Qq, and that ft = at f. Then by the definition of p\ 
as the density of the image measure of Q\ one obtains that the right hand side 
of (|39|) is bounded by 

^ [ l ge Pl(q)\ft(q)\dqdt 

Using the bounds (|19~D (with p t replaced by p\) and (|34|) (which holds on C^ e ) we 
finally obtain that 



P({? e : \Di(q)\ > K}) < 



i f/ I 

K Wo in 2 



9|Vt(g)p 



dt 



dqdt + pj^ J^\V^ t (q)\ 2 dq dt 



(40) 



The bracket on the r.h.s. is (|35|). By the Assumptions A3 and A4, (|35|) is finite 
and hence the r.h.s. of ( (4~0|) goes to zero uniformly in e as K — > 00. Thus we have 
established Lemma 13. (A □ 
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4 Bohmian mechanics and self-adjointness 



4.1. In this subsection we shall discuss the necessity of certain assumptions 
under which we have established global existence of the Bohmian particle mo- 



tion (cf. Theorem [3.1| and Corollary |3.2|) . We shall investigate in particular the 
assumptions concerning self-adjointness of the Hamiltonian. 

By Corollary |37| we obtain global existence if the Hamiltonian is the form sum 
H + V, and if the potential V satisfies certain conditions leading in particular 
to the Hamiltonian's being bounded from below. These conditions on the Hamil- 
tonian guarantee in particular that Assumption A4 of Theorem |3.1| is satisfied. 
In the case of one particle moving on the half line Q = (0, oo), we shall prove, 
without invoking A4, global existence for a certain class of potentials for arbitrary 
self-adjoint extensions, which furthermore may be unbounded below. 

Theorem 4.1 Let VI = (0, oo), H = L 2 (VL), and suppose V G C°°(Q) is such 
that Hq + V is in the limit point case at infinity (see for example $3^/)- Let H be 
an arbitrary self-adjoint extension of (H + V)\cg°(n); and let ip G C°°(H) with 
\M\ = 1. Then P(r < oo) = 0. 



It follows for example from Theorem X.8 in [32] that if V(r) > — kr 2 for r > c 
with c, k > 0, then H + V is in the limit point case at infinity. 

Consider as an example the potential V(q) = —c/q 2 with c > large enough: 
The Hamiltonian H = H + V is in the limit circle case at 0, in the limit point 
case at infinity, and unbounded above and below (cf. for example [j32"|). Thus 
by Weyl's limit point-limit circle criterion there is a one-parameter family of 
(similarly unbounded) self-adjoint extensions of H |cg°(Q) f° r a h of which, by Theo- 
rem |4.1| , Bohmian mechanics exists uniquely and globally for P-almost all initial 
values. 

The proof employs a new definition of the particle dynamics in one dimension 
which extends the solution to @ and is interesting in its own right. (In fact, this 
definition extends the Bohm motion, defined by ([[]) and @, to an equivariant 
motion for all ip G L 2 \) Let Qt{qo) be defined implicitly by 

i-Qt{qo) ri° , 



Mq)\ 2 dq= / \Mq)\ 2 dq 

-oo J — oo 

Qt(qo) is well-defined if 

F(q,t) := T \^ t \ 2 dx 
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is strictly monotonic in g. This is the case except at extended intervals with 
ipt = 0, where F(-,t) has a plateau. To define Qt{%) globally for q Q G IR, set for 
example 

Q t (q ) := min{g : F(q, t) = F{q , 0)} (41) 
(and Q t (q ) = -oo if F{q Q , 0) = 0, Q t (q ) = oo if F(q , 0) = 1). 



Proof. From Lemma [TT] we obtain that ip G C°°(f2 x IR). Therefore, using the 
continuity of the scalar product and the L 2 -differentiability of 1 1 — > -0*, 



F(q,t) = \ip t \ dx = (%>,,] V>t, A) 

(where (•, •) denotes the scalar product in 7i = L 2 (Q)) is jointly continuous and 
differentiable. Clearly F{0,t) = 0, ]jm q ^ OB F{q,t) = 1, and dF/dq = |^ t (g)| 2 . 
Moreover, 

dF(q,t) f«d\ilH 



at 



I -^tf~ dx = + liir k 3t{c) = -Jt(q)- 

Jo at c^o 



Here the existence of lim c ^o jt{c) follows for ip G C°°(H) from partial inte- 
gration of f c (ip*(Hip) — (Hifj*)ip) dx and Schwarz's inequality; the value for 
lim c ^ jt{c) — follows from the symmetry of H together with the fact that 
Hindoo jt{d) = 0, which holds because H is in the limit point case at infinity. 
(See for example [|37|j.) 



For all t and all q G Go = Q\jV , let Qt(<lo) be defined by (S). It follows from 
the implicit function theorem that 1 1— > Qt{qo) is continuous and differentiable for 
(q ,t) such that MQM) ± 0, with dQ t /dt = j(Q t )/\MQt)\ 2 = v**(Q t ), i.e., 
Qt solves the differential equation on Q — (O x IR) \ Af. It remains to show 
that for P-almost all initial go, r (<?o) = sup{s > : Qt(qo) G G for all t < s} is 
infinite, i.e., fl2|) has global solutions for almost all initial values. 

Now it is obvious from this definition that Qt{qo) G Q for all t and all go G Go- 
(Qt(<lo) — corresponds to F(g ,0) = 0, Qt(lo) = oo to F(g ,0) = 1, and for 
go £ Go, F(q o ,0) G (0,1).) Moreover, by the L 2 -continuity of t — > tp t , we have 
that for < T < oo and g G Go, inf <t<T Qt(lo) > and sup < t < T Q t (qo) < 
oo, i.e., the trajectories cannot run into the (only) possible singularity of the 
potential S = {0} or to infinity in finite time. Thus it remains only to control 



the probability of hitting M , for which Lemma |3.5| does the job. We omit the 



details. □ 
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4.2. One might now wonder whether we have global existence of Bohmian me- 
chanics for any self-adjoint Schrodinger Hamiltonian (without assuming A4). 
This is quite trivially wrong, as is easily seen by considering free motion on 
the interval Q = (0, 1). There are self-adjoint extensions of #o|cg°(n) with 
j(0) = 7^ 0- (Similarly one might consider potentials on Q = (0, oo) such 
that Hq + V is in the limit circle case at infinity. ) This corresponds to an incoming 
flow at 0, balanced by an outgoing flow at 1 (or the other way round) so that 
the total probability is conserved (a situation which can of course be identified 
with a motion on a circle). Typically, the particle will reach the boundary of Q, 
so that almost sure global existence in the sense of solutions of the differential 
equation (Q) fails. However, the motion is quite trivially extendible in such a way 
that the trajectories are piecewise solutions of the differential equation: when the 
boundary of Q is reached they jump to the other end of Q. \ip\ 2 then remains an 
equivariant measure. This motion can be described by replacing ( fill) by 

Q t (q ) :=mm{q:F(q,t) = F(q ,0)} 

with 

F(q,t) = (F(q,t)- £j a (0)ds^ (mod 1) 

[Another possibility to define a global motion in this case is to use the unmodified 
(PI). This provides then an example of a deterministic dynamics completely 
different from (and not an extension of) the Bohmian dynamics, (([2j) is replaced 
by the nonlocal form dQ/dt = — jt(0))/|V'| 2 ) for which, however, |-?/>| 2 remains 
equivariant. With this motion, particles do not jump from 1 to or the other 
way round. (However, they might all run through nodes!)] 

In fact, we expect generally that self-adjointness guarantees (possibly discontin- 
uous) extendibility of the Bohmian motion in such a way that |^| 2 is an equivari- 
ant measure. This is suggested by the fact that the symmetry of the Hamiltonian 
leads to 

lim [ / (fUq) ■ u)ds + / (j**(q) ■ u)ds] = 0, 

using integration by parts (Green's identity) 

/ ip*(Hip)dq- [ {H^*)^dq = -ih [ (f-u)ds, 

JM JM JdM 



27 



for M = JC n \S s . The vanishing of the integrals over the absolute flux yields global 
existence of Bohmian mechanics: In finite time the singularities and infinity are 
not reached. The flux balance from self-adjointness alone suggests extendibility 
of the motion: Some parts of the singularities (or infinity) may act as sources, 
others as sinks. 



4.3. For a wider perspective on this matter let us consider a Schrodinger Hamil- 
tonian H on a domain where it is not (essentially) self-adjoint, i.e., where the 
boundary conditions are too few or too weak. Then, first of all, the time evolu- 
tion of wave functions is not unique: There are infinitely many different unitary 
evolutions (corresponding to the different self-adjoint extensions), and there are 
also semi-groups for which \\ipt\\ is not conserved. The (essential) self-adjointness 
of H is equivalent to Ker(H* ± i) = {0}, so that if H is considered on a domain 
where it is symmetric but not self-adjoint, then H* has imaginary eigenvalues. 
Together with the (space) regularity for eigenstates of the elliptic operator H* 
(assuming sufficient regularity for the potential V) we thus obtain classical solu- 
tions of Schrodinger's equation with exponentially decreasing or increasing norm. 
Since p = \ip\ 2 still holds on X (cf. the paragraph around Equation (fT^D), those 
solutions lead with positive probability to catastrophic events. 

This possibility is not that far-fetched: The Hamiltonian for one particle in a 
Coulomb field V(r) = — 1/r considered on the "natural" domain Cg°(IR 3 \ {0}) 
is not essentially self-adjoint and hence the time evolution of the wave function 
is not uniquely defined |22|, |19[ . There are many properties that mathematically 
distinguish the self-adjoint extension usually regarded as "the Coulomb Hamilto- 
nian" from other possible extensions. However, we do not know of any convincing 
(a priori) physical argument for "the Coulomb Hamiltonian" unless one accepts, 
for example, that the Coulomb potential is a "small perturbation" of the free 



Hamiltonian ||20|| , or that "in reality the singularity is smeared out." Of course, 
if we require that Bohmian mechanics be globally existing, then, as we have ar- 
gued above, only self-adjoint extensions are possible. But among all self-adjoint 
extensions Bohmian mechanics seems not to discriminate: While our Corollary 
^T^ applies only to the form sum (which is "the Coulomb Hamiltonian"), it is 
heuristically rather clear (or at least plausible) that Bohmian mechanics should 



exist globally and uniquely for all the other self-adjoint extensions of H 



C-(R 3 \{0}) 
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as well.[] 

Nonetheless, discussions about the "right" (unitary or contractive) evolution, 
i.e., about the "right" boundary conditions, as for example in the case of strongly 
singular potentials like the 1/r 2 potential (see [ll], do gain now firm 

ground by taking into account the actual behavior of the particles: Whether 
or not we should consider the Bohmian particle to be caught at the origin is a 
matter of the physics we wish to describe: whether or not particles disappear in 
the nucleus. An axiom, or dogma, of self-adjointness of the Hamiltonian (or equi- 
valently of unitarity of the wave function evolution) appears quite inappropriate 
from a Bohmian perspective — even though the importance of self-adjointness is 
profoundly illuminated by this perspective! 

Moreover, the particle picture of Bohmian mechanics naturally yields an in- 
terpretation of the current j as a current of particles moving in accordance with 
the density \ip\ 2 - In this way, boundary conditions for self-adjointness of the 
form j = at the singularities or j(in) = j(out) may be viewed as "arising 
from Bohmian mechanics." For example, the outcome of a detailed analysis 
of self-adjoint extensions of H on the half line (0, oo) — there is a one parame- 
ter family of self-adjoint extensions Hq, the respective domains being defined 
by ^'(O) / ^(0) = a, a real, or ip(0) = (a = oo) — is easily guessed from the 
point of view of Bohmian mechanics by demanding that either v^(0) = 0, i.e., 
Im(V>'(0)/V>(0)) = 0, or that |^(0)| 2 = 0. 



4.4. We wish to conclude with some remarks on the general Hilbert space de- 
scription of orthodox quantum theory viewed from the perspective of Bohmian 
mechanics. We have discussed the fact that Bohmian mechanics is well defined, 
i.e., trajectories exist globally and uniquely, for typical initial values and for 
wave functions which are C°°-vectors of the self-adjoint Hamiltonian H. The set 



9 The singularity of the radial current at may be estimated for if; = J^i m fim(r)Yi m (0 , 4>) S 
0^ o D(H r j) ® Ki, where, for angular momentum I, H r> i is the radial part of H and K\ is the 
corresponding eigenspace of the angular part of — A, as follows: H r ifi m £ L 2 (TR + ,r 2 dr) 
implies for the worst behavior of fi m as r — > that fi m ~ r a with a > 1/2 for I > 1 resp. 
a = — 1 for I = 0. Therefore the worst behavior of the the radial current jf = ^Im as 
r — * is jf <~ r~( 3 / 2 ) +e (using the fact that the radial current at vanishes on V(H r _o) for all 
self-adjoint extensions, cf. the proof of Theorem 4J). Thus we should have that JU r \j^-u\ds = 
f$2 Mt \r 2 duj ^ r( 1 / 2 ) +e — > as r — > 0. For a proof of global existence along the lines of Theorem 
|3.l|, it is necessary also to control the time change of the radial current. However, the global 



existence for the one-dimensional problem (Theorem 4.1) suggests that this should be possible. 
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C°°(H) is dense and invariant; however, it is most likely not a residual in the 
norm topology of the Hilbert space L 2 (Q), i.e., it is presumably not a "generic" 
set, and it furthermore depends on the Hamiltonian. One might now wonder how 
Bohmian mechanics can be taken as the basis for the quantum formalism (as has 
been claimed — see [[11]]) if the former cannot even be defined for a really "fat" 
set of wave functions. And since, as we have seen, Bohmian mechanics yields a 
natural understanding of the (spirit of the) meaning of the self-adjointness of a 
Schrodinger Hamiltonian, the question should be even more puzzling. The an- 
swer is, of course, that the embedding of Bohmian mechanics into a Hilbert space 
structure is a natural but purely mathematical device. Indeed this answer is (of 
course, in disguise) commonly accepted — though maybe not as loudly stated: No 
physicist believes that a generic L 2 -wave function (in the residual sense) results 
as the "collapsed" wave function from a preparation procedure. The state space 
of physical wave functions ip is not the Hilbert space 7i = L 2 (Q) but more or 
less the space of classical, smooth solutions of Schrodinger's equation, for the 
analysis of which the L 2 -norm and hence the Hilbert space structure is of critical 
importance. 

Other aspects of this embedding are commonly taken more seriously: for exam- 
ple, that observables are self-adjoint operators on 7i. While we do not wish here 



to enter into a general discussion of this question (see flllj), we would like once 
again to comment on the self-adjointness of the Hamiltonian H. The importance 
of this property is certainly not that "measured energy values must be real" but 
lies rather in Stone's theorem: H acts as the generator of a one-parameter unitary 
group U t , which gives the time evolution of states ip t = U t ip (or of observables 
A t = Ut l AU t ), and hence must be self-adjoint by Stone's theorem. Why should 
the time evolution be unitary? Simply because the norm \\ipt\\ must be invariant, 
so that the total probability is conserved. 

We conclude with some remarks about effective descriptions. We first note that 
restrictions of configuration space such as described in the last paragraph of Sec- 
tion 4.3 (with a freedom in the boundary condition) are perhaps best understood 
physically as arising as a limit of a sequence of (moderately realistic) potentials 
V£ tending to "V = for q > 0, V — oo for q < 0" in a suitable way — such that 
Hq + V® — > Hq in an appropriate sense. This problem is analyzed in [[3^, [I], ||], 
and in || the convergence of the Bohmian trajectories in this limit is derived. 

In other physically interesting but complex situations we may have an effective 
description involving a Hamiltonian which is self-adjoint but not of Schrodinger- 
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type ,0 so that the probability current j may fail to be of the usual form [cf. 
Eq. (H)], or where there may in fact be no local conservation law at all for the 
probability density j^j 2 . 

For an example with nonstandard current j, consider the self-adjoint shift 
operator H c = —ihcV, where c is a constant with the dimension of a velocity. 
ipt(q) = e~ ltHc / h ipo(q) = ipo(q — ct) describes translation without "spreading." 
This Hamiltonian may perhaps arise in a limit in which the spreading of the 
wave function, induced by the Laplacian, can be neglected. In any case, the 
corresponding current is j c = c\ip\ 2 , and the obvious candidate for the "Bohm 
motion" in this case is v = j c /p = c, not flU). 

It is conceivable that an approximation procedure leading to an effective Hamil- 
tonian like H c , when applied to Bohmian mechanics, also converges to a deter- 
ministic limit. If so, then = c would be the natural guess for the motion in 
this limit. 

For a Hamiltonian with no local conservation law for probability there is of 
course no "Bohm motion" generalizing ([!]). 
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6 Appendix: On the regularity of ip 

Lemma 6.1 Assume A1-A3, and let if)t{o) = c~ ltH ^ h ipQ{q). Then there exists a 
function ip G C°°(Q x IR) such that for allt e IR ip(q, t) = ifitio) f or almost all q. 
(ip is a classical solution of Schrodinger's equation.) 

This fact is presumably folklore knowledge to experts in PDE's, but since 
we could find no suitable reference — and since it does not appear to be well 

10 The modeling of physical situations leads often to idealizations which are very singular. In 
Newtonian mechanics one considers for example singular evolutions induced by "hard walls" 
confining a particle or by elastic collision between hard spheres. 
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known among mathematical physicists — we shall supply a proof. (Hunziker |]T8| 
has established space-time regularity of ip for potentials which are bounded, have 
bounded derivatives, and are C°° on IR d for ip in Schwartz space. Also, regularity 
(in space) of eigenf unctions (for sufficiently regular potentials) is well known |[32| .) 

Proof. We apply standard methods of elliptic regularity (see, for example, |33| ) 
to the elliptic operator L on Q x JR 

r) 2 N ft 2 B 2 

From A3, ip t e C°°(H), and therefore the functions n)t := H n ip t (= e~ itH/h (j) nfi ) 
are in L 2 (Q) for all n and t. With this definition, formally 

14 = 02 + 0i, and L0 n = n+2 + </> n+ i (42) 

To apply the theorem of elliptic regularity, we need to show a) that ip and 4> n 
are locally L 2 in Q x H, therefore locally in the Sobolev space W° (we refer to 
the definitions and theorems of [Q, where however W n is written H n ) in Q x 1R, 
and b) that (021) is satisfied in the distributional sense on Q x 1R. Then, by 



repeated use of Theorem 8.12 in |3l| we obtain that ip (and <p n ) are locally in W n 
for all even (positive) integers n. Then by Sobolev's lemma i/j is indeed (almost 
everywhere equal to) a C°°-function on Q x 1R. (The space-time set of measure 
on which ip has to be corrected indeed splits into t-slices that are of measure 
for all t. This is a consequence of L 2 -continuity of t > ip t .) 

a) This is an easy consequence of Fubini's theorem if ijj and 4> n are jointly 
measurable in q, t. ip t is measurable in q 

(t/j t e L 2 (n)) and the map t i—* ipt 
resp. t 1 — > </>„, ; i is weakly measurable (indeed much more is true, namely strong 
differentiability). Then by a theorem of Bochner and von Neumann || joint 
measurability of ipt{q) and <p n ,t{q) in follows in the following sense: There 
exist functions ijj and <j) n which are jointly measurable in q, t, and for all t ip(q, t) = 
ipt(q) and <f) n (q, t) = n ,t(<z) f° r almost all q E Q. In the following, we shall denote 
il>(q, t) and <p n (q, t) by ip(q, t) and 4> n (q, t) or again by ipt(q) and n>i , as convenient. 

b) First one convinces oneself that ip and <p n satisfy Schrodinger's equation in 
the distributional sense, i.e., for all test functions / G C^°(fl x H), 

d . 

dt. 



-i% J (-f)iP dqdt = J (Hf)tP dq 
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This follows by looking at the function G : 1R — > IR, t i— > / f(q,t)if)(q,t)dq = 
(fti ^t)) where (•, •) denotes the scalar product in L 2 (Q). G has compact support, 
and its derivative is seen to be 

dG(t) 1 ,„ rr , , ,df* , . 

by the continuity of the scalar product (•, •) and the weak (in L 2 (Q)) differentia- 
bility of if) t . The same holds with if> replaced by <p n for all n. Since furthermore 
ft G T^{H) for all i, and therefore (f*,Hifj t ) = (Hf£,if) t ), if) and n are indeed 
weak solutions of Schrodinger's equation (and the self-adjoint operator H on 
T>(H) agrees with the operator on distributions defined by H), i.e., we've arrived 
at 

ih— = Hip = 0i and ih—^ = H(p n = (p n+1 
at at 

weakly, and therefore ( [42|) indeed holds in the distributional sense (on Q x 1R) . □ 
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